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We develop an efficient and general method for optimizing the microwave field that achieves 
magnetization switching with a smaller static field. This method is based on optimal control and 
renders an exact solution for the 3D microwave field that triggers the switching of a nanomagnet 
with a given anisotropy and in an oblique static field. Applying this technique to the particular 
case of uniaxial anisotropy, we show that the optimal microwave field, that achieves switching with 
minimal absorbed energy, is modulated both in frequency and in magnitude. Its role is to drive 
the magnetization from the metastable equilibrium position towards the saddle point and then 
damping induces the relaxation to the stable equilibrium position. For the pumping to be efficient, 
the microwave field frequency must match at the early stage of the switching process the proper 
precession frequency of the magnetization, which depends on the magnitude and direction of the 
static field. 

We investigate the effect of the static field (in amplitude and direction) and of damping on the 
characteristics of the microwave field. We have computed the switching curves in the presence of 
the optimal microwave field. The results are in qualitative agreement with /i-SQUID experiments 
on isolated nanoclusters. The strong dependence of the microwave field and that of the switching 
curve on the damping parameter may be useful in probing damping in various nanoclusters. 

PACS numbers: 75.50.Tt,75.75.-n,75.10.Hk 
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I. INTRODUCTION 



Fine magnetic clusters offer tremendous challenges 
both in the area of fundamental science and practical 
applications. The main reasons for this impetus are the 
novel features related with their small size, such as the 
possibility of high density storage, short-time switching^ 
and fast read-write processes. On the other hand, the 
small size is a drawback in many regards. The energy 
barrier in these systems is too small to ensure a reason- 
able stability, in a given energy minimum, that is neces- 
sary for practical applications at room temperature, e.g., 
magnetic recording. This is the problem of superparam- 
agnetism. A possible way out would be to use materials 
with high anisotropy and thus ensuring a high energy 
barrier. A consequence of this is that high values of the 
writing (or switching) fields are required. However, it 
is still unclear how to devise such high fields operating 
on the scale of nanoclusters while avoiding the ensuing 
noise. In order to keep the size small, the energy barrier 
high, and the switching field small, other routes are ex- 
plored and a promising one among them is provided by 
microwaves. Microwave-assisted magnetization switch- 
ing in various magnetic systems, such as thin films, has 
been investigated by many groups^,. In fact, we have 
at hand a more general and fundamental issue, namely 
the problem of getting a system out of an energy mini- 
mum by nonlinear resonance. This has previously been 
addressed in many areas of physics and chemistry, espe- 
cially in the context of atomic physics. For example, Liu 
et al4 have studied the dissociation of diatomic molecules 
by a chirped infrared laser pulse and showed that this 
process requires a much lower threshold laser intensity 
to achieve dissociation. The quantum regime^ has been 



studied in terms of energy-ladder climbing and gives very 
similar results. Experimental evidence of this process has 
been provided by the dissociation of HF molecules using 
a sub-nanosecond frequency modulated laser pulsed. 

According to the classical theory of auto-resonance or 
the quantum theory of ladder-climbingS"^^, exciting an 
oscillatory nonlinear system to high energies is possible 
by a weak chirped frequency excitation. Moreover, trap- 
ping into resonance followed by a (continuing and stable) 
phase-locking with the drive is possible if the driving fre- 
quency chirp rate is small enough. It has also been shown 
that a slow passage through and capture into resonance 
yields efficient control of the energy of the driven sys- 
tem. Incidentally, an important theoretical result is that 
the precise form of the time dependence of the oscillating 
field is not essential for the process to succeed. 

For nanoclusters, it has been shown in previous works 
how a monochromatic microwave (MW) pulse can, by 
means of a non-linear resonance, substantially reduce the 
required static field needed to reverse the magnetization 
of an individual nanoparticle. Indeed, it has been demon- 
strated using the /U-SQUID technique on a 20-nm cobalt 
particle^ that adding an MW field with given rising time, 
duration, and frequency, on top of the static (DC) mag- 
netic field, the switching of the nanocluster is possible 
at a static field lower than the Stoner-Wohlfarth (SW) 
switching field and within a time interval of the order 
of a nanosecond. The switching curves, or the so-called 
SW astroids, obtained in these measurements present 
some irregular features implying that the reduction of the 
switching field is not uniform, as is the case with thermal 
effects. The global features depend on several physical 
parameters, such as the MW field pulse duration, its ris- 
ing time, and its frequency, the DC field amplitude, and 



2 



the damping parameter—. The SW astroid obtained with 
its peculiar features seems to bear the fingerprints of the 
nanocluster and its underlying characteristics such as its 
potential energy. Moreover, the strong dependence of 
these features on the damping parameter might be used 
to estimate the latter in such clusters. 

On the theory side, several works have been devoted to 
the understanding of the magnetization dynamics, and 
in particular its reversal, under the effect of a time- 
dependent magnetic field. The theoretical work may be 
divided into two kinds. The first deals with the effect of 
a given MW field with a given polarizatio n^°i^^ while the 
second seeks optimal strategies for achieving the magne- 
tization switching^. In particular, a few works, e.^.— , 
assume a given dynamics for the magnetization and at- 
tempt to determine the MW field that realizes it. 

In the present work, we use a general method borrowed 
from the optimal control theory^^^^^ and apply it to the 
switching of a nanocluster. This method renders an exact 
solution for the MW field vector necessary for the switch- 
ing of a nanomagnet with a given potential energy (com- 
prising anisotropy and an oblique static field) . The stan- 
dard formulation of this method consists in minimizing 
a cost functional using the conjugate gradient technique. 
The latter is known to be a local-convergence method 
and thus renders a solution that is rather sensitive to 
the initial guess. In order to acquire global convergence 
and thus have solutions for the MW field that are sta- 
ble with respect to a change in the initial conditions, we 
have supplemented the conjugate gradient routine by a 
global search using the Metropolis algorithm and simu- 
lated annealing. Then, we have applied our algorithm to 
a nanomagnet in the macrospin approximation with uni- 
axial anisotropy and oblique DC magnetic field. We have 
investigated the effect of the latter (both in direction and 
magnitude) and of damping on the characteristics of the 
MW field. Then, we computed the limit-of-metastability 
(or switching) curves for different (small) values of damp- 
ing. 



II. METHOD OF OPTIMAL CONTROL 
APPLIED TO NANOMAGNETS 

One of our objectives here is to develop a general 
method that allows us to solve the following inverse prob- 
lem: what is the optimal time-dependent magnetic field 
under which the magnetization of a nanoparticle, with 
given potential energy, switches from a given initial state 
to a given final prescribed target state? After formulat- 
ing this method we apply it to the case of a macrospin, 
or a nanoparticle in the SW approximation, in an energy 
potential composed of uniaxial anisotropy and a Zeeman 
contribution from an oblique DC field. 

The method we propose is borrowed from the opti- 
mal control theory. The main idea is to start with an 
arbitrary MW magnetic field hAc (t) with its three com- 
ponents /iaC' = 2;, y, z, that we call the control field, in 



addition to the static magnetic field and anisotropy field. 
We then determine hj^c{t) that triggers the switching 
of the cluster's magnetization between two given states 
within a prescribed interval of time. 



A. Model 

Consider a nanomagnet in the macrospin approxima- 
tion where its magnetic state is represented by a macro- 
scopic magnetic moment m — ii^s, where /i^, is its mag- 
nitude and s its direction with |s| = I. In this ap- 
proximation, the relevant terms entering the energy E of 
the nanomagnet are the magneto-crystalline anisotropy 
and the Zeeman energy. The applied DC (static) field 
Hdc is assumed to point in an arbitrary direction = 
Hdc/^dc- Using the convention /Iq 1 so that the 
magnetic fields are expressed in Tesla, one then defines 
the effective field 
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and writes the damped Landau-Lifshitz equation in the 
Gilbert form (LLE) that governs the dynamics of s, as- 
suming that the module fj.^ remains constant. 



I ds 

jdi 



= — s X Hcff — a s X (s X Hoff ) 



(2) 



where 7 ~ 1.76 x 10^^ (T.s)^^ is the gyromagnetic factor 
and a the phenomenological damping parameter (taken 
here in the weak regime). 

We measure all applied fields in terms of the anisotropy 
field 



Ha - 



2KV 



(3) 



and in particular we define the reduced effective field 

1 ss 

hcff = —Hcff = (4) 

Ha OS 

with £ = E/ {2KV). In terms of hcff LLE becomes 
ds 

— -s X hcff - a s X (s X hcff) (5) 
dr 

where r = t/tg is the dimensionless time and tg — 
l/{'jHa) the characteristic scaling time of the system. 
For instance, for a cobalt particle of 3 nm diameter— 
with K ~ 2.2x 10^ J.m'^, ~ 3.8 x lO^^o ^.m^ we 
have Ha ~ 0.3 T and ~ 1.9 x 10"" s. 



B. Formulation of the optimal control problem 

a. General procedure The idea here is to introduce 
a control field hAc(''') = Hac/ Ha and then seek its op- 
timal form that allows for driving the magnetic moment 
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direction s from the given initial state s^*^ at time = 
into the desired final state s*^-''-' at the given observation 
time r / . Accordingly, we replace in the LLE ^ the (de- 
terministic) field hoff by the total (time-dependent) field 



C(r) 



loff ■ 



lAc(,Tj 



(6) 



This results in the following equation of motion, which 
will be henceforth referred to as the driven LLE (DLLE) 



-s X C{t) — as X (s X C(''')) ■ 



(7) 



The field hAc(''') is then determined through the mini- 
mization of a cost functional which, in the present case, 
may be written as 



J'[s(T),hAc(T)] = X s(t/)-s 



(8) 

The first term measures the degree at which the magnetic 
moment switching is achieved and vanishes in the case 
of full switching. The second term is quadratic in the 
driving field and is thus proportional to the absorbed 
energy. The parameter 77, called the control parameter, 
allows us to balance the second condition with respect to 
the first. 

Therefore, the problem of optimal control boils down 
to minimizing the cost functional ([5]) along the trajectory 
given by DLLE ([7]). More explicitly, this amounts to 
solving the following problem 



in^^[s,hAc]-i s(r/)-s'^^ " + f /drhic(r) 



s = — s X ^ — as X (s X ^) , Te[0,T/] 



s(0) = sW. 

(9) 

An optimal solution of this problem is characterized 
by the first order optimality condition in the form of the 
Pontryagin minimum principle (PMP)^-. These condi- 
tions are more conveniently formulated with the help of 
a Hamilton function which may be in the present case 
written in the following form 



H[s(T),A(r),hAc(r)] 



(10) 



A(r) • {— s X ^ — as X (s X ^)} , 



where A(r), called the adjoint state variable [see be- 
low], is a Lagrange parameter introduced to implement 
the constraint and thereby render s(t) independent of 
hAC (''')• The PMP then states that solving the problem 
^ is equivalent to solving the following boundary prob- 
lem (i.e., the Hamilton- Jacobi equations with boundary 



conditions) 

fs=M, s(0) = sW, re[0,r;], 
^ = -^: A(r/)=s(r;)-s(/), (11) 
= 

ShAC 

The last condition is also equivalent to the vanishing of 
the gradient of the cost functional J- in Eq. (|8]). It yields 
the equation 



sn 

ShAC 



= rj Hag + sxA — asx (sxA), (12) 



which is used to compute the variation in the cost func- 
tional, that is 



5F= dr 
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(13) 



In general, this problem is highly nonlinear and con- 
sidering, on top of that, the non-linearity of the Landau- 
Lifshitz equation, it is not possible to find analytical solu- 
tions. Consequently, we resort to numerical approaches. 
The advantage of this formulation is manifold: i) the 
MW field hAc(''') is obtained in 31?, i.e., one obtains 
the three functions of time /i^(-.(r),a = x,y,z; and for 
any potential energy (anisotropy, DC field, etc), ii) the 
final time tj, and the absorbed power (second term in 
Eq. ([8|)) can be adjusted; the latter may be achieved by 
tuning the control parameter 77, iii) one can generalize 
this treatment to many-spin problems^^ and also include 
thermal effects. 

One of the most efficient techniques for (numerically) 
solving such a minimization problem is the conjugate- 
gradient method. However, the drawback of this method 
is that it is a local-convergence method, which means that 
the solution it renders is strongly dependent on the initial 
guess. We overcome this inconvenience by supplement- 
ing the method by a global search using the Metropolis 
algorithm with random increments and then proceed by 
the technique of simulated annealing. 

For numerical calculations, we have discretized the 
boundary- value problem (jlip by subdividing the time in- 
terval [Ti = 0,T/] into N time slices 



Tn = Ti + n X At, 



0, , 



,iV-l, 



where 



At = 



N -1 ' 



Then, using the notation v„ — v(t„) for a vector v, 
Eqs. ([71 [51 [T^ and the equation for A, become 
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-1 = S„ + Ar X [-S„ X Cn -aSnX (s„ X Cn)] 1 s(Ti) = S 



(^) 



^=2 



2 

2 / y ~AC,n' 
ji=0 



A„_i = A„ - At X A„, \f = sn^i - s^f\ 



At X [77 hAc,n + s„ X A„ - a s„ X (s„ x A„ 



(14a) 

(14b) 

(14c) 
(14d) 



The explicit expression for A„ in (fT4b) depends on 
the energy potential [see below for the case of uniaxial 
anisotropy] . 

We may summarize the numerical procedure as follows, 
i) for a given initial guess of the control field hAc(Oi 
first solve the state equation (IT4k ) forward in time using 
the initial condition, and then evaluate the cost func- 
tional (|14b ). ii) the solution obtained for s is then used 
for the backward (since the condition now is at tf) inte- 
gration of the equation ([Mb) for A, iii) with the trajec- 
tories of s and A thus obtained we compute the gradient 
(I14d ). The numerical subroutines are standard and can 
be found in Ref. [l^- We emphasize that obtaining the 
control field amounts to solving for 3 x iV variables. 

h. Uniaxial anisotropy In the case of uniaxial 
anisotropy with oblique static field the energy of the 
nanomagnet reads (in units of the anisotropy energy 
2KV) 



£ = -hnc (e/i • s) - i(s • n)^ 



(15) 



with K and n being the anisotropy constant and easy 
axis, V the nanomagnet volume and hnc = Hue /Ha- 
The effective field explicitly reads [see Eq. |4] 



hoff = /iDC e,, + (s • n) n. 



(16) 



From the second equation in pip we obtain the explicit 
equation for A 

A = C X -'^ + a [C X (A X s) + A X (C X s)] 

+ [A- (s X n + a s X (s X n))] n (17) 

and in Eq. (jl4b ) we now have 

■'^ri = Cn X A,i 

+ a [Cn X (A„ X 

Sn) + A„ X (Cn X S„)] 

+ [A„ + a (A„ X s„)] • (s„ X n) n. 



III. RESULTS 

In the present work, we have considered the case of a 
nanomagnet with uniaxial anisotropy and oblique static 
field. Unless otherwise stated, the latter is applied in the 
yz plane making an angle of 170° with respect to the easy 
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Figure 1: Optimized MW field (upper panel) and the corre- 
sponding spin trajectories (lower panel). The inset is a 3-D 
plot of the spin trajectory on the unit sphere. 



axis (z axis). Its reduced magnitude is hue — 0.5, corre- 
sponding to a field magnitude Hue — 150 mT. The initial 
position and target states s^'^ and s'-^\ which correspond 
respectively to the metastable equilibrium state and sta- 
ble equilibrium state, are computed numerically. The ob- 
servation time is t/ = 600 (i.e., tf ~ 11.4ns). The damp- 
ing parameter is a = 0.05 and the control parameter 77 
has been set to 0.01. The static field, damping parameter 
and observation time have been varied and their effects 
studied [see later on] . For simplicity, we have taken a lin- 
early polarized MW field, i.e., hAc(^) = h\c{i)^x- This 
choice also suits the experimental setup^. 

In Fig. [1] we have plotted the optimized MW field 
magnitude Hxc (t) = Ha h^c (0 where t is the time in 
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Figure 2: Instantaneous frequency of the optimized and fil- 
tered MW field of Fig. [1] The other parameters are the same 
as in Fig. \T\ 



seconds, together with the components of the magnetic 
moment, i.e., Sq {t) ,a = x,y,z. First, we note that the 
amplitude of the MW field is rather small as it does not 
exceed 15 mT, which is 10 times smaller than the static 
field. Moreover, the summed magnitudes of the DC and 
MW field are smaller than the SW switching field for the 
chosen DC field direction (about 200 mT). This shows 
that, in the presence of a MW field, magnetic switching 
is achieved at a smaller DC field. Second, the striking 
feature is that the MW field is modulated both in am- 
plitude and frequency. Its frequency is a slowly varying 
function of time in the stage that precedes switching, as 
can be seen in Fig. [5] Third, as is hinted to by the dashed 
vertical lines, the extrema in the MW field and the spin 
components Sy{t) and Sz{t) match at all times before 
switching. This simply implies that the magnetic mo- 
ment is phase-locked to the MW field. All these features 
agree with the predictions of the classical auto-resonance 
or the ladder-climbing quantum theory, as summarized 
in the introduction. 

The instantaneous frequency has been obtained after 
passing the optimized MW field through the Butterworth 
filter and then applying the Hilbert transformation2ii2^. 
As can be seen, for short times the instantaneous fre- 
quency oscillates around the approximate value /o ~ 
4.1 GHz. This initial frequency is simply the FMR fre- 
quency given by 



/] 



FMR — 



2tt 
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where h^^l n = hcgi^*^ 



s(^) = /iDC (e„ • s(»)) + (s(*) • n)2 
is the effective field ([T6| evaluated at and then pro- 
jected onto the initial position s*^'). As the magnetic mo- 
ment approaches the saddle point the frequency decreases 
rapidly and eventually vanishes when the magnetic mo- 
ment crosses the saddle point into the more stable energy 
minimum. 



In Fig. 1 it is seen that the time span comprises three 
stages (for the set of physical parameters considered): 
1) Nucleation stage (up to 5.4ns). The MW field re- 
mains almost zero and the magnetic moment remains in 
the metastable state. 2) Driven precession (from 5.4 ns 
to 9.7 ns). Here the MW field and the magnetic moment 
are synchronized. At each procession cycle, the MW field 
hooks up the magnetic moment and pushes it upwards 
in the energy potential towards the saddle point. This 
is the phase-locking process mentioned in the introduc- 
tion and observed above. This is indeed possible because 
the frequency chirp rate is small as can be seen in Fig. 
[2] for 5.4 ns < t < 9.7 ns. The MW field thus compen- 
sates for the effect of damping that tends to pull the 
magnetic moment back towards its initial position. At 
around 9.7 ns, the magnetic moment crosses the saddle 
point. 3) Free relaxation: from 9.7 ns onward, the mag- 
nitude of the MW field dwindles and the synchronization 
with the magnetic moment is lost. We note that at the 
saddle point the precession reverses from being counter- 
clockwise to clockwise as the magnetic moment switches 
to the lower half sphere. 

Numerical tests show that the MW field can be re- 
placed by zero during the nucleation and free relaxation 
stages without noticeably affecting the trajectory of the 
magnetic moment. This implies that the most relevant 
part of the signal is that during the driven precession; the 
role of the MW field is thus to drive the magnetic mo- 
ment towards the saddle point. Next, the damping takes 
up to lead it to the more stable energy minimum. Dur- 
ing the driven precession the frequency of the MW field 
and the precession frequency of the magnetic moment are 
similar. Consequently, the magnetic moment switching 
can be viewed as a resonant process: the pumping by 
the MW field is efficient when its frequency matches the 
frequency of the magnetization (phase- locking). 

The same calculation has been carried out with the 
same sampling time but different values for the total ob- 
servation time tf. The results are shown in Fig. [3] If the 
total time is larger than an effective time of 6 ns, sim- 
ilar values are obtained for the cost functional and the 
curves /iac(^) can be matched after a time shift. As was 
discussed earlier, this effective time corresponds to the 
sum of the time of driven precession and that of free re- 
laxation. This result implies that the nucleation stage 
can be suppressed without affecting the final optimized 
MW field. However, if the total time is too short, the final 
value found for the cost functional is higher (i.e. not fully 
minimized). Indeed, we see in Fig. [3] (uppermost panel) 
that the stage of driven precession is shortened and the 
shape of the control field changes so as to achieve a faster 
switching and thus comply with the switching-time con- 
straint [first term in Eq. (jS])]. 

The effect of varying the amplitude of the static field on 
the MW field is shown in Fig. |4l We see that the shape of 
the MW field envelop remains the same, apart from the 
fact that the smaller the static field, the more symmet- 
rical is the MW field. This shows that for a higher field 
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Figure 3: Optimized MW field obtained with four different 
total times tf. 
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Figure 4: Optimized MW field obtained for different magni- 
tudes of the static field hoc, in the same direction making an 
angle 6 = 170° with respect to the anisotropy easy axis. 



/iDC) the energy potential is less symmetrical. Moreover, 
as /iDC is increased the energy barrier is lowered and the 
MW field required to achieve switching is smaller. Again, 
the initial frequency of the oscillations matches the FMR 
frequency of the system; when /idc increases, the latter 
decreases. The cost functional was found to be propor- 
tional to the energy barrier between the saddle point and 
the metastable minimum. Hence, when the energy bar- 
rier is higher, more energy has to be injected in order to 
overcome it. The same study has been carried out upon 
varying the direction of the static field. 

We have also investigated the effect of varying the 
damping parameter a on the MW field. The results are 
summarized in Fig. [5l We see that the intensity of the 
field increases with a, which is compatible with what was 
suggested earlier, namely that the role of the MW field is 
to compensate for the damping effect. This effect is sim- 
ilar to what happens with a rubber band; the more you 
stretch it the harder it becomes to do so. Moreover, the 
effective duration of the MW field, which mainly corre- 
sponds to the driven precession period, decreases when a 
increases. We note that, on the contrary, the initial fre- 
quency of the oscillations is independent of a [see Fig. [S]- 
This result can be understood qualitatively if we suppose 
that, at any time, the MW field exactly compensates for 
the effect of damping. The spin dynamics is then gov- 
erned by the undamped LLE and the magnetic moment 



precesses with its proper frequency, which is independent 
of the damping parameter. At short times, since the pre- 
cession angle is small, this precession frequency is equal 
to the FMR frequency. 

As discussed in the introduction, one of the objectives 
of investigating the magnetization switching assisted by 
MWs is to achieve an optimal switching with smaller DC 
magnetic fields than it would be necessary without MWs. 
This means that applying the DC field in a given di- 
rection and varying its magnitude one determines the 
switching field (or the field at the limit of metastability) 
at which the magnetization is reversed. This is the SW 
astroid. Due to the energy brought into the system by 
MWs, the field required for switching is smaller. This has 
been nicely demonstrated using the /i- SQUID technique 
on a 20-nm cobalt particle^. The most striking feature 
of the SW astroid obtained by these measurement is its 
jaggedness. In other words, the reduction of the switch- 
ing field is not uniform and presents a kind of "fractal" 
character. The global features depend on several physi- 
cal parameters, such as the MW field pulse duration, its 
rising time, its frequency, the DC field amplitude, and 
the damping parameter. In the present work, and in the 
particular case considered here, namely that of uniaxial 
anisotropy, we first wanted to check whether this reduc- 
tion of the switching field is recovered by our optimal- 
control method. Furthermore, we address the question 
as to whether the SW astroid may be used as a finger- 
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Figure 5: Optimized MW field obtained for different values 
of the damping parameter a. 
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Figure 6: Instantaneous frequency of the MW field optimized 
for several values of the damping parameter a. Dotted line: 
FMR frequency. 



print of a given nanocluster. More precisely, the ques- 
tion is vifhether a given SW astroid can provide us with 
specific information about the corresponding cluster, like 
its energy potential and the physical parameters such as 
damping. 

Accordingly, we check whether an MW field h^(^(i), 
which is optimized in the presence of a reference ap- 
plied DC field h.^Q{t) with given direction and magni- 
tude, e.g. hue — 0.5 and an angle of f70° with respect 
to the easy axis, can still induce magnetization switching 
in the presence of another DC field, with different di- 
rection and/or magnitude. To answer this question, the 
MW field h^c(i), was used in the driven LLE © and the 
calculation of the switching field was performed for sev- 
eral intensities and directions of the static field h leading 
to the switching curves in presence of h^^ {t) as shown in 
Fig. [71 As can be seen, the magnetization switching oc- 
curs only inside the golf-club-shaped green area [see Fig. 
[7] (left)]. In the black area, the pumping by the MW field 
is inefficient and switching does not occur. This curve is 
in agreement with the experimental data of Ref. |^. 

The shape of the green pattern can be explained based 
on qualitative arguments about the frequency and mag- 
nitude of the MW field. As has been seen previously, in 
order to achieve the switching, the MW field must fulfill 
the following conditions: i) it must be synchronized with 
the proper precession frequency of the magnetization; so 
at short times its frequency must match the FMR fre- 



quency of the system, and ii) the injected energy, must 
be sufficient to overcome the energy barrier between the 
metastable minimum and the saddle point. 

For any magnitude or direction of the field Hdci both 
the FMR frequency and the energy barrier can be com- 
puted numerically [see Fig. [7] (right)]. In the black area 
the value of the FMR frequency is the same as for h.^^. 
In the hatched area the energy barrier is lower than for 
the hpQ. Outside the black zone, the MW field is not 
synchronized with the precession frequency of the sys- 
tem: the switching can not occur. Outside the hatched 
area the injected energy is not sufficient to overcome 
the energy barrier. Consequently, the switching is only 
achieved in the intersection between both areas. Indeed, 
comparing with Fig. [7] (left), this intersection matches 
more or less the green zone, where the switching occcurs. 

Next, we optimize the MW field h.^^{t) for the ref- 
erence DC field hpQ(t) with magnitude hue = 0.5 and 
angle of 170° with respect to the easy axis, and damping 
a'^ = 0.05; then we compute the SW astroid for other 
values of a, in the presence of the same DC and MW 
fields. The results are shown in Fig. [5] 

We see that the shape of the switching area strongly 
depends on the damping parameter a. The largest green 
area is found for a = a^. Then, as a increases the green 
area shrinks and vanishes for a > 0.12. Indeed, for high 
values of a, the MW field is not strong enough to com- 
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Figure 7: (a) Switching curve computed in the presence of 
the MW field hAc(^)- The red cross indicates the amplitude 
and direction of the DC field for which the MW field was 
optimized. The area in green is where switching has been 
achieved, the black area is where there is no switching, and 
in the white area the static field is higher than the switching 
field (i.e. beyond the metastability region), (b) In the black 
area the FMR frequency is the same as for the reference DC 
field hnc • In the hatched area the energy barrier between the 
metastable minimum and the saddle point (computed numer- 
ically) is smaller than for the reference DC field hpc- 

pensate for the effect of damping. The same phenomenon 
is observed for small values of a, in which case the MW 
field "overcompensates" for the effect of damping and 
thereby the energy can not be pumped into the system 
in an efficient manner. 



IV. CONCLUSIONS AND OUTLOOK 

We have developed a general and efficient method for 
determining the characteristics (pulse shape, duration, 
intensity, and frequency) of the MW field that triggers 
the switching of a nanomagnet in an oblique static mag- 
netic field. We have applied the method to the case of 



uniaxial anisotropy and investigated the effect of the DC 
field and damping on the optimized MW field. We have 
shown that our method does recover the switching field 
curves as observed on cobalt nanoclusters. It remains 
though to investigate the origin of the "fractal" charac- 
ter observed in the measured switching curves. 

We have shown that the MW field that triggers the 
magnetization switching, while minimizing the absorbed 
energy, can be efficiently calculated using the optimal 
control theory. According to our results, the optimal MW 
field is modulated both in frequency and in magnitude. 
The role of this MW field is to drive the magnetization 
towards the saddle point, then damping leads the mag- 
netic moment to the stable equilibrium position. For the 
pumping to be efficient, the MW field frequency must 
match the proper precession frequency of the magnetiza- 
tion, which depends on the magnitude and the direction 
of the static field. Moreover, the intensity depends on the 
damping parameter. This result could be used to probe 
the damping parameter in experimental nanoparticles. 

The present method is quite versatile and can be ex- 
tended to other anisotropics. It could also be used to 
study the dynamics of nanoclusters in the many-spin 
approacbiS. In this case one will probably have to 
deal with a nonuniform MW field, especially if surface 
anisotropy is taken into account^^. One may then study 
switching via internal spin wave excitations and the effect 
of the MW field on the corresponding relaxation rate^— . 
Thermal effects can also be accounted for by adding a 
Langevin field on top of the DC and MW fields. In this 
case, it will be interesting to investigate the interplay be- 
tween the MW field and the Langevin field and to figure 
out when these two fields play concomitant roles. 
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